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We propose a phenomenological theory for topological crystalline insulators with time reversal 
and Ca symmetries. First, we introduce a fictitious space and transformation of electromagnetic 
field operators. This transformation leaves the speed of light unchanged but changes the elementary 
charge to \ple. Then we formulate the theory of topological crystalline insulators in terms of 
transformed fields in this fictitious space as 3D BF theory containing 7r-flux excitations. It is known 
that a 3D BF theory with half flux quantum excitations describes low energy properties of time 
reversal invariant insulators. By making an inverse transform we recover the effective field theory 
in original space. It turns out that this field theory contains quarter flux quantum excitations. 

I. INTRODUCTION 

In a recent paper—, Fu showed that there are 3D topological band insulators characterized by a Z2 index protected 
by time reversal and C4 symmetries. These insulators have gapless surface states, 2D massless fermions with quadratic 
spectrum, only on one of their surfaces (the one which is perpendicular to the C4 symmetry axis) . Other surfaces of 
these insulators turn out to be gapped. 

Three-dimensional time reversal invariant insulators^— are characterized by the following topological response*^ 
the axion term (c = 1) 

9e 2 
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with 9 = (2n — l)ir in the non-trivial phase. Since the axion term ([T]) is related to a surface quantum Hall effect 
with the conductivity ^ %- , and since the Hall conductivity of massless 2D fermions with quadratic spectrum is twice 
as bigger than that of massless 2D fermions with linear spectrum, one may expect that electromagnetic response of 
topological crystalline insulators of RefQ] is given by an axion term with 9 = 2tt. We show that this is indeed the case 
by employing the technique for computation of the axion term described in the paper @. 

There is a topological field theory description of 3D time reversal invariant insulators^ which is given by 3D BF 
theory 

Cbf = ^-e^ pX a,8 v b pX + L £ ^ Advb 
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where k = 2 . This theory can explain all of the main properties of 3D time- reversal invariant insulators: Z2-ness, 
axion term, gapless edge theory. One may wonder is there an analogous theory for crystalline insulators protected by 
point group symmetries X 

In constructing such a theory one is faced with some difficulties which are clear already at this stage of discussion. 
One of them is related to the fact that C4 symmetry of the underlying microscopic model can not enter the topological 
field theory in an explicit way and should be hidden. Second is related to the fact that an axion term with 9 = 2tt 
is equivalent to an axion term with 9 = due to 2-7T periodicity^&^iiS which makes it difficult to explain ^-ness of 
topological crystalline insulators in the context of topological field theory. 

In this paper, we construct a theory that allows to circumvent difficulties mentioned above. This theory explains 
Z2 classification of topological crystalline insulators not in terms of original fields, but in terms of fields transformed 
to a fictitious space. It turns out that unit charge in this space is different from that of an electron charge and as a 
consequence it follows that 2-7T axion electrodynamics transforms to 7r axion electrodynamics. This is crudely how Z2 
classification is explained. This theory also contains traces of C4 symmetry because it supports quarter flux quantum 
excitations. 

II. COMPUTATION OF THE TOPOLOGICAL RESPONSE 

Consider the following hamiltonian 

H(k x , k y , k z ) = h(k x ,k y ) ® t z + M ■ 1 (g) r x - k z ■ 1 ® r y 
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where 

k 2 1 . 1 

h(k x , k y ) = — + ^(k 2 x - k 2 y ) ■ a z + —k x k y ■ a x (2) 

If M > for z > then (in the r space) z-dependence of the wavefunction is 




Analogously, if M < for z < then 




This means that there is a 2D masslcss chiral fcrmion with quadratic spectrum ([2j localized near the surface which 
is perpendicular to the symmetry axis and which separates two models with different sign of mass M. 

There are two symmetries of the hamiltonian, time reversal operation T which coincides with complex conjugation 
K, and C4 rotation in the (k x , k y ) plane U = a y 1, (UT) 2 = —1: 

T H(^k X: ky : k z *)T = U( k Xl k y: k z ^) : U U(k x , k yi k z ^)U = H(k y , k X: k z *) 

We will compute the axion term assuming that h(k x ,k y ) = ■ 2 -^ (k 2 — k 2 ) ■ a z + -^-k x k y ■ cr x , where m\ = rri2, using 
the approach given in the appendix of Ref. @. Calculation can be easily generalized for unequal masses mi 7^ m%, but 
unfortunately we could not do the calculation for generic h, Eq. ©. Lagrangian in the euclidean space-time is 

^(7^ - M)ip 

where 

70 = l®f x , 1l=(T z ®T yi "j2=<Jx®T y , 73 = l(g)T z , 75 = a y (8 T y = 70717273 



d a = d t +iA a , d 1 = —(k 2 ~kl), d 2 = - — {k x k y + k y k x ), d 3 = d z +iA z , /,-,..„ <>,.„ I iA x>y 

M — > Me* 75 interpolates between hamiltonians with different sign of mass when a varies from to 71", violating 
time reversal symmetry for intermediate values of a. If a slowly varies in space such that there is a domain wall on 
which the mass M flips the sign, then the action on one of the sides of the domain wall will contain an additional 
term with respect to the action on the opposite side. This additional term is an anomaly due to the change of the 
Jacobian under the transformation ip' = e lQ75 / 2 V>- This anomaly can be computed by a proper regularization, using 
Fujikawa method, and is given by 

m \wf ^n( k \M) 2 /M 2 \k)}} (3) 

where the function /(s) satisfies the requirements 

/(0) = 1, /(oo)=0, s/'(s) = when s = 0,oo 

It is easy to show that (7^) 2 = d^d^ + 1/4^, 7„][d M , d v ]. Since tr{7 5 [7 M , 7„][7 P , = 16e^ pcr , in the absence of 
the gauge fields, = 0, one obtains (k\ — (d^d^)\k) = k 2 + k 2 + (k 2 + k 2 ) 2 /Am 2 . To obtain a non-zero result in 
© in the limit M -> 00 we need to expand f up to the second order in ["y^ , d„]/M 2 and retain the terms in 

[d^.d^Wdp.da] which contain the operators d 2 and d 2 . Namely, the terms that should be retained are 

[d , d 3 ] = iE z , [di,d 2 ] = -^B z (d 2 + d 2 ) + ... (4) 

111 A u 



[d Q ,d X ] = —{Eydy-E X d X ) + 



[d 3 ,d 2 ] = -{B y d v -B x d x ) 



(5) 
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[d ,d 2 ] = ^ i (E x dy + E y d x ) + ..., [d 3 ,d 1 ] = -±-{B y d x + B x d y ) + ... (6) 

Then (J3|) becomes 

m j ^• 4 -^W{[^,di][rf2,d3]-Mo,d 2 ]Mi,d 3 ] + [4,d3]Mi,d2]}|fc>r( s ) = *^ J ^■4- J ^-(Bvy(kl+kl)f"( s ) 

where s = k 2 . + k\ + (k x + k 2 ) 2 /4m 2 . Calculating the integral, which does not depend on the specific choice of the 
function /, we get the axion term i^^BE) with 8 = 2n. 

III. TRANSFORMATION OF ELECTROMAGNETIC FIELD 

In this section we introduce a fictitious space (denoted in the following by ~) and a transformation of electromagnetic 
field which is motivated by the previous discussion, Eqs.((3HB]). The aim is to reformulate the theory in this fictitious 
space and in terms of transformed electromagnetic field operators. 

It is more convenient to work in momentum space. Momentum operators in ~ space arc given by 

K x = — ^r— ^, K V = ^L, K z = k z (k y >0) (7) 

and with reverse sign for k y < 0. From this definition one can see that the absolute value of the momentum does not 
change since K 2 + K 2 = k"]_ : and that the vector K covers " momentum space twice. This means that electromagnetic 
field operators in ~ space split into two fields (generally, every field that can be transformed into ~ space will split 
into two fields). Wc will divide kj_ plane into two regions: 1) sgnfcj, ~ sgnk y (I = 1) and 2) sgnk x = — sgnfc. y (I = 2). 
The relation between k and K integrations is given by 



1=1,2 

Scalar field in ~ space is defined according to relations 



2 J d 3 k= J d 3 K (8) 



W(K) = M, (k x ,k y >0); ^)(K) = ffi, >0) (9) 



and (/>"'(— K) is defined as complex conjugate of Eqs.©, such that in accordance with © we have 

J ^(k)0(-k)d 3 fc = J W (K)0 (/) (-K)d 3 /i 



1=1,2 ' 



Transformation formulas for vector potential must be consistent with Eq. (O. They are defined as follows 

AW(K)= kxAx{k) ~ kvAy{k \ iW(K)= ^ Aa(k) + fca ^ (k) , i(0(K) = ^M (10) 
V2kj_ V2fcj_ v2 

and A^(— K) as complex conjugate of (|TU)) . One can see from Eq. (fTTJj) that when vector potential in original space 
is a pure gauge, then it is pure gauge also in ~ space, which is quite reasonable. Transformation formulas for electric 
and magnetic fields can be obtained by substituting (O and ([TUf into equations Ei = —dtAi — di<j> ,Hi = EijkdjAk- 
The result is that E and E, H and H are related to each other exactly by the same formulas as A and A, Eq. (|10p . 
and E")(— K), HW(-K) by their complex conjugates: 

^) (K) = fe^(k)-^(k) ^( 0(K) = ^(k) + ^(k) ^) (K ) = ^ k l (11) 
V2k ± y V2k ± 2 V2 

^ WO-WO ^ )(K) = ^(k) + W) ff(o (K) = ^ (12) 

V2fcj_ y V2fcj_ v2 



4 



It is also possible to incorporate source terms (current j^) into the theory, transformation properties of which are 
analogous to ([5]). ([ID]). 

One of the consequences of Eqs. (fTTj) and (IT21) is that scalar products in two spaces differ by a factor of 2. For 
instance, E(K)E(— K) = E(k)E(— k)/2, and analogous formulas are true also for the scalar product of the pairs H(K) 
and H(-K), E(K) and H(-K). This leads to equalities 

f E 2 d 3 x = E / [E (0 ] 2 rf 3 x, f H 2 d 3 x =J2f [H (i) ] 2 d 3 x 

J 1=1,2^ J 1=1,2 

which mean, in particular, that lagrangian density of electromagnetic field does not change its form and, therefore, A 
is a true electromagnetic potential in ~ space. Electric and magnetic fields (| 1 1 1) and ([12")) defined in this way satisfy 
Maxwell's equations with the same speed of light c. 

Coulomb gauge condition kA(k) = transforms to KA(K) = 0. So, it can be checked by direct computation 

that the commutation relations^ [Ei(k), Aj (— k')] = 4ni (5ij — ^pr 2 -^ • (27r) 3 <5(k — k'), which are satisfied in quantum 

theory in the Coulomb gauge, are transformed to [.e{ (K), Af\-K')] = Am (s tj - ■ (2tt) 3 5(K - K.')6 U >. This 

means that transformation formulas (| 10[) arc correctly defined also in quantum theory. 

To see how the axion term transforms one needs to consider 9 varying in space. If 9 = 9(z) then, rewriting the 
integral in momentum space, it can be shown that 



J #(z)EH d 3 x = Y,J ^)E ( °H (0 d 3 a 



So one can write 

EH = L E W H« (13) 



M , JM 



IV. SOME OTHER PROPERTIES OF THE FICTITIOUS ELECTRODYNAMICS 

One of the consequences of Eq. © is that the scalar potential of a static point charge q, 0(k) = Airq/k 2 , becomes 
(f)(K) = 4:Trq/K 2 , q = q/V2, i.e. scalar potential of a static point charge q/V2- More generally, a simple correspondence 
can be made between electric field distributions in the two spaces when the field configuration has axial symmetry. 
However, this does not mean that elementary charge in ~ space e is given by 1/V2 units of electronic charge. 

In the same way, one can obtain the transformation q m — > q~ m = q m /\/2, <p — > (f> — for a magnetic pole q m 

and a magnetic flux vertex <j) parallel to z axis. 

There is one more configuration which has simple transformation properties, a skyrmion. Consider the simplest 
skyrmion in 2+ ID 

4>o{r,ip) = cos/(r), <j> x {r, ip) = sin/(r) cos<^, <j> y {r, ip) = sin/(r) sin^ (14) 
This can be written in more convenient form as 

</> = cos/(r), (t> x = d x smf(r')dr', (f> y = d y I sin f{r')dr' . 

Jo Jo 

Transformation formulas for <pi are analogous to (fT0|) (wc drop the factor ^ because it does not affect the final 
result), and in coordinate space we obtain 

4> (r, tp) = cos /(?'), 4>x\r, ip) = sin /(r) cos ifi® (r, <p) = sin f(r) sin tp, 1 = 1,2, (15) 

where r and (p are polar coordinates in ~ space. 

When a skyrmion is coupled to a doublet of two-component fermions through the interaction lagrangian ga z (<pf), 
it acquires quantum numbers, charge and spinJ^— For the skyrmion coupled to Dirac fermions induced fermion 
current is given by the topological current of the skyrmion (we put number of flavours Nj = 1) 



Jfj, — —-£ tiv\£ijknid v njd\nk 



(16) 
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where m = fa/\(p\, and spin is determined by the Hopf term 

S Hop f = ™ I d 3 x,J^C^ (17) 



where the gauge field is related to the current through the relation J M = e^xd^Cx- Charge is given by the 
temporal component of (|16p 



Q = I d 2 x^-Eij k nid x njdyn k 



1 

The skyrmion configuration (|14j) has charge Q = e and spin s — ^. Wc note that spin can not be treated according 
to quasiclassical equations (fT0|) . Correct quasiclassical treatment of the spin of a skyrmion is as in 

Rcfslilr{l3 

However, we are interested in skyrmions in fermion models with quadratic band touching. These were studied 
ir>i£ and the main result can be formulated as the doubling of quantum numbers compared to fermions with linear 
spectrum. Thus, we have that the skyrmion (|14[) coupled to fermions with quadratic spectrum has charge Q = 2e 
and spin s = 1. This is analogous to doubling of the 9 term that we have found in Sec. [TTJ After transformation to 

space, the Hopf term Suopf = 2iir J drxJ^C^ of such a skyrmion will be equally distributed between skyrmions 

(US): S Hop f = Snlpf + Sjffopf' ^Hopf = i7T I d^xJ^C^, where J M = ^s^ v \Sij k hid v njd\hk, hi = fa/\4>\ (we note that 
fermionic charge current is = \f2eJ lJ ). This means that skyrmions in ~ space will have half of the spin of the 
skyrmion in original space. In the case under consideration, their spin is ^, i.e. they are fermions. 

We saw that the correspondence Suopf = Shov} ~^^Hopf holds for the skyrmions (fT4"|) and (fT5|) . but it can be shown 
that it is more general. One can see that fa depends continuously on fa. If we continuously deform fa then Suopf 
does not change, since it is topological invariant. Analogously, S^ op f does not change under continuous deformations 

of fa. So SHopf = Sff op f + SHopf remams valid for continuous deformations of the skyrmions ([T4]) and ([15]) . 

Thus, we have established that transformation of spin-1 bosonic skyrmions with electric charge 2e to ~ space results 
in two different types of spin-i fermions each having charge y/2e. Suppose that we interpret this as e = \/2e. Since 
9 term is proportional to the square of elementary charge, this means that the value of the 0-anglc in ~ space is two 
times smaller than in original space. This idea is pursued further in the next section from different point of view. 

Now we consider an example for illustration, an axially symmetric charge distribution p(r) rotating around the 
symmetry axis z with small angular velocity uj (oja <C c, a is characteristic size in transverse direction). We have 
j x = —uid y g(r), j y = u>d x g(r), j z = 0, where g(r) = f£ p(r')dr' . Transformation to ~ space gives p^ l \r) = -^p(f), 

ji 1 ^ = —u!d y g(r), jy = uid x g{r), I = 1,2 where g(r) = J Q r p(r')dr' . This means that angular velocity in " space is the 
same as in original space. In particular, rotation around z axis in original space of an axially symmetric configuration 
results in rotation by the same angle in ~ space. Wave function of a bosonic skyrmion (S z = 1) acquires phase factor 2ir 
after 2ir rotation in original space. However, in " space we have phase factor it for each sector I, such that 2tt = n + ir. 
Together with the fact that spinor (S z = 1/2) acquires phase factor tt after 2tt rotation, this means that the skyrmion 
(|15[) in ~ space is a fermion. 



V. BF THEORY 



In this section, after briefly reviewing 3D BF theory, we turn to the study of the following possibility that a 
system in ~ space is described by a 3D BF theory with k = 2. We will show that transformation of electromagnetic 
field operators defined in section Mil is consistent with quantization of charge in ~ space in y/2 units of elementary 
charge. This in turn means that a 9 = 2n axion electrodynamics becomes a 9 = tt axion electrodynamics in ~ space. 
Transformation of this model to the original space gives a BF theory with k = 4 and we argue that this is an effective 
field theory for topological crystalline insulators with C4 symmetry. 

A. Review of 3D BF theory 

3D BF theory contains two statistical fields a M and coupled to quasiparticle and vortex currents 

Cbf = ^-e^ x a^d v b p x+fa^ + Js^ft^ 

47T Z 
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The first term imposes non-trivial statistics between quasiparticlcs and vortices: when a quasiparticle j is rotated 

k 



around a vortex E the full wavefunction of the system changes by a phase ^ £^ In the absence of quasiparticles and 
vortices lagrangian can be written as 



k f 

Lbf = ^- / (PxcLibi + a (dibi) + b i(eijkdjak) 

where 6; = \sijh^jk- Integrating out ao and &oi gives the constraints dibi = and Eijkdjdk = which can be resolved 
on a torus by taking^ 

en = y~ + diA, bi = ^- + £ijkdj(k 
where A and C arc periodic functions, LiSi is the volume of the system. The final lagrangian contains only a, and 6.; 

Lbf = -^-aidtbi 

Z7T 

It is possible to construct gauge invariant quantities from cii and bi, Wilson loops and Wilson surfaces^ 

A t = e le& \ Bi = e i<l>oli 

This means that a,; = a,i+ 2 ^- and bi = b~i + ^ which is a manifestation of quasiparticle and vortex number quantization. 



B. Effective field theory description of topological crystalline insulators 

a and A appear in the BF theory in the combination a + A, so a transforms in the same way as the electromagnetic 
potential. However, a on a compact surface should be treated differently because it is not possible to define a as the 
K — > limit of the continuous expression. Instead we write 

l-W^E^f, k = \ (18) 
l 

where = af^ = cLi/j, &j = b^ = jbi for some numeric constant 7. Since Wilson loops and Wilson surfaces in 
space are given by 



A ■ — p iga i R — pi<t>abi 



e<p = ecpa 



it follows that e = ye. Thus bi = j 2 ■ bi/y, which means that transformation formulas for b are given by that of a 
with an extra factor of 7 2 , such that 

7 2 M a^ = l^/.afW (19) 



2tt J m ' 2ir — JM 



For the two Eas.(fT8f and (|19[) to be consistent with each other 7 must be equal to This means that e = V2e. 
The difference between transformation of a and b can be viewed as the consequence of charge-flux duality, since 
transformation formulas for a (b) differ from that of momentum K by additional factor l/v2 (y/2). 

The full Lagrangian in ~ space together with time-reversal symmetry breaking surface terms is given by the expres- 
sion (fc = 2) 

E ^ vpX ^dJ>% + ^ vpX AV~d v b% + ^P%aPd p Af> (20) 

Lagrangian (|20p is two copies of the model that was proposed in the Ref. HI as an effective field theory for 3D topological 
insulators. Transformation of the model (f2"Uf to the original space gives a model with the Lagrangian density (fc = 4) 

Ctci = -e^a^bpx + -e^ px A p d u b pX + ^ £ ^ x d p a u d p A x (21) 

7T 7T 47T 
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An interesting consequence of Eq. (|21[) is that vortex excitations in this model carry 7r/2-flux, i.e. quarter flux quantum. 
This can be associated with C4 symmetry, since due to C4 symmetry of the microscopic model it is expected that 
quarter flux quantum must play a special role in the theory 

Transformation of A, a and b under time-reversal is the same as in original space: (Aq, — A4), (So, — 5j) and {—bo, bi) 
after time reversal, so the action (|2ip is time-reversal invariant. C4 rotation in original space interchanges the sectors 
with I = 1 and / = 2 with each other. 

If one integrates out a and b from Eq. (|2ip then the remaining term is a 2tt axion term. Eq. (|13p . together with 
the fact that unit charge in ~ space is e = \/2e, gives that 

e 2 9 f l 2 9 ^ f ~ ~ -6 

— / EH = — V / EH, 8=- 

2nhJ M 2 

Thus a 9 — 2-7T axion term transforms into two copies of 9 — n axion terms. These two copies arc different, because 
they involve Fourier components of electromagnetic field with different momenta. So one can not combine these two 
BF theories into a single BF theory with 9 = 2tt. These two theories should be treated as independent unless there 
are other terms mixing them with each other. 

Z2 classification of topological crystalline insulators can be explained as it is done for 3D topological insulators, by 
applying the arguments of RefQ to the model ([2U| . Situation described here is, in some sense, similar to the calculation 
of Zi invariants of topological insulators, where Brillouin zone is divided into two halfs, related to each other by time- 
reversal symmetry^ The difference is that here we deal with an effective field theory instead of microscopic model: 
by a suitable transformation, topological field theory is divided into two and a logic, familiar from the study of 3D 
topological insulators^ is applied to each half. 

Eq . ((2"Tj) is consistent with general consideration of effective field theories of topological insulators.— It was shown 
inl23l that before specifying the properties of space-time, the lagrangian density of topological insulator is of BF type 
with two unknown coefficients 

A F (e^ px a^b pX + e^A^bpx) + C F e^"%a v d p A x 

This theory assumes that collective excitations in the bulk of a topological insulator is given by neutral fcrmions, 
which couple to external fields by electric and magnetic dipole moments. It is reasonable to assume that collective 
excitations in topological crystalline insulators also can be described by this method. 



VI. DISCUSSION 



We have presented some arguments that the theory (|2II) might be an effective field theory for recently proposed 
topological crystalline insulators^ in the case of C4 symmetry. However the discussion has not been rigorous and there 
are some points that should be cleared. For example, transformation of momentum operators Q was defined only in 
the continuous limit and we do not know what are the corresponding relations for a compact manifold. Such relations 
are necessary for a rigorous theory. We also did not explain why the model in " space is given by two copies of the 
model that was proposed as a topological field theory for time-reversal invariant topological insulators^ One of the 
reasons for this choice was that this theory can easily explain Zi property and we did not know any other theories. 
Though it was not a priori stated that the theory must contain excitations with quarter flux quantum, we obtained 
that there are such excitations making quite reasonable propositions about the structure such a theory should possess. 
This last circumstance makes these propositions more plausible. 

The other question is related to gapless surface excitations. It was shown in the ReflU how one can account for 
gapless surface excitations of 3D topological insulators in the context of BF theory. Electronic degrees of freedom on 
the surface are constructed using 2+1D bosonization scheme proposed in the Ref l24H25l In principle, some analogous 
procedure can be done in the case of fermions with gapless quadratic spectrum in the continuous limit, by reducing 
them first to fcrmions with linear spectrum. However it is more proper to study fermions with quadratic spectrum 
from the beginning in view of the following circumstance. In 2+1D in the continuum limit, internal energy of a 
single two-component fermion with gapless quadratic excitations is equal to internal energy of a boson with gapless 
quadratic excitations due to the identity 




k 2 d 2 k f k 2 d 2 k 



analogous to equivalence of internal energies of bosons and fermions in 1+1 Di2& No such simple relation exists for 
2+1D fcrmions and bosons with linear spectrum and bosonization of Dirac fermions with linear spectrum in 2+1D, 
as it was described i n 24 ' 25 , is not fundamental^ 
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